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Abstract 


Geometry  optimizations  and  normal-mode  analyses  of  three  conformers  of 
1,3,5-trinitro-s-triazine  (RDX)  are  performed  using  second-order  Moller-Plesset  (MP2)  and 
nonlocal  density  functional  theory  (DFT)  methods.  The  density  function  used  in  this  study  is 
B3LYP.  The  three  conformers  of  RDX  are  distinguished  mainly  by  the  arrangement  of  the  nitro 
groups  relative  to  the  ring  atoms  of  the  RDX  molecule.  N02  groups  arranged  in  either  pseudo- 
equatorial  or  axial  positions  are  denoted  with  (E)  or  (A),  respectively.  The  axial-axial-equitorial 
(AAE)  conformer  has  Cs  symmetry  and  is  the  structure  in  the  room-temperature-stable  crystal 
(a-RDX).  The  axial-axial-axial  (AAA)  and  equitorial-equitorial-equitorial  (EEE)  conformers 
have  C3v  symmetry,  a  symmetry  consistent  with  vapor  and  (5-solid  infrared  (IR)  spectra.  The 
AAE  and  AAA  conformers  are  studied  at  die  MP2/6-31G*,  B3LYP/6-31G*,  and 
B3LYP/6-31 1+G**  levels,  and  the  EEE  conformer  is  studied  using  the  B3LYP  density  function 
and  the  6-3 1G*  and  6-31  1-Kj**  basis  sets.  The  geometric  parameters  and  IR  spectra  of  the  AAA 
conformer  are  in  good  agreement  with  experimental  gas-phase  and  P-solid  data,  supporting  the 
hypotheses  derived  from  experiment  that  the  AAA  structure  is  the  most  probable  conformer  in 
vapor-phase  and  P-solid  RDX.  The  B3LYP/  6-31 1-fG**  structures  and  simulated  IR  spectra  are 
in  closest  agreement  with  experimental  data.  The  MP2/6-31G*  structures  and  spectra  are  in 
poorest  agreement  with  experiment. 
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1.  Introduction 


Advances  in  the  development  of  nonlocal  density  functional  theory  (DFT)  [1-5]  and  in  computer 
architectures  have  allowed  for  reliable  electronic-structure  investigations  of  large  polyatomic 
molecules.  We  have  subjected  the  well-studied  explosive,  1,3,5-trinitro-s-triazine,  commonly  known 
as  RDX,  to  both  DFT  and  ab  initio  treatments.  The  results  presented  here  provide  atomic-level 
information  about  RDX  conformers,  as  well  as  indicate  the  suitability  of  current  theoretical 
treatments  to  systems  such  as  these.  In  this  work,  we  determine  the  geometries  of  three  conformers 
of  RDX  and  characterize  them  through  normal-mode  analyses  using  electronic-structure  methods. 
Comparisons  of  structural  parameters,  vibrational  frequencies,  and  simulated  infrared  (IR)  spectra 
against  measured  properties  are  given. 

There  have  been  attempts  at  treating  RDX  with  electronic-structure  and  semi-empirical  theories 
[6-7];  however,  the  highest  level  of  theory  used  for  geometry  optimizations  reported  in  these  studies 
is  the  SCF-MO  (Hartree-Fock  [HF])  level  using  basis  sets  ranging  from  STO-3G  to  4-21G  [6-7]. 
These  levels  of  calculations  can  provide  a  crude  approximation  to  optimized  structures;  however, 
known  deficiencies  in  the  theory  beg  for  further  theoretical  treatments  to  provide  a  reliable  prediction 
of  the  system.  An  example  resulting  from  the  deficiencies  in  the  HF  theory  is  seen  in  the  SCF/4-21G 
calculations  of  Coffin  et  al.  [7],  in  which  SCF  geometry  optimizations  were  attempted  for  four 
conformers  of  RDX.  The  four  conformers  of  RDX  differ  mainly  in  the  position  of  the  nitro  groups 
relative  to  the  ring  atoms.  The  ring  atoms  are  arranged  in  the  chair  conformation.  The  conformers 
were  labeled  according  to  axial  (A)  or  pseudo-equatorial  positioning  (E)  of  the  nitro  groups  about 
the  ring.  We  adopt  the  same  nomenclature  for  the  conformers  in  this  work.  All  nitro  groups  of  the 
axial-axial-axial  (AAA)  conformer  occupy  axial  positions,  and  all  nitro  groups  of  the  equitorial- 
equitorial-equitorial  (EEE)  conformer  occupy  pseudo-equatorial  positions.  Both  of  these  conformers 
belong  to  the  C3v  point  group.  Two  nitro  groups  occupy  axial  positions  on  the  axial-axial-equitorial 
(AAE)  conformer  and  the  remaining  nitro  group  is  in  the  pseudo-equatorial  position.  For  the  AEE 
conformer,  two  nitro  groups  occupy  pseudo-equatorial  positions,  and  the  remaining  nitro  group  is 
axial  [7].  The  AAE  and  AEE  conformers  belong  to  the  C5  point  group.  Figure  1  illustrates  all  of 
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Figure  1.  Structures  of  the  AAE,  AAA,  and  EEE  Conformers  of  RDX.  Atom  Labels  on  the  Two-Dimensional  Projection  of  the 
RDX  Molecule  Are  Consistent  With  the  Internal  Coordinates  in  Tables  1  and  2.  For  Clarity,  the  Hydrogen  and  Oxygen 
Atoms  Are  Not  Dlustrated  in  the  Three-Dimensional  Representations. 


these  conformers,  except  for  the  AEE  species.  For  clarity,  we  have  not  shown  the  hydrogen  atoms 
or  the  oxygen  atoms  of  the  nitro  groups  for  the  AAA,  EEE,  and  AAE  conformers  in  this  figure. 
Normal-mode  analyses  for  each  of  the  SCF/4-21G  optimized  structures  indicated  that  only  the  AAA 
conformer  had  all  real  vibrational  frequencies  and  was  therefore  the  only  stable  species  predicted  at 
this  level  of  theory.  The  remaining  conformers,  each  of  which  had  at  least  one  equatorial  nitro 
group,  had  imaginary  frequencies,  indicating  they  were  not  stable  structures.  Neutron  diffraction 
measurements  of  solid  a-RDX  (the  form  that  is  stable  at  room  temperature)  provided  the  crystal 
structure  and  atomic  arrangements  of  the  molecules  in  the  crystal  [8].  The  molecular  geometry  of 
a-RDX  is  consistent  with  the  AAE  conformer.  This  indicates  that  either  the  SCF/4-21G  level  of 
theory  is  not  sufficient  to  correctly  describe  this  conformer,  or  that  the  crystal  field  stabilizes  the 
AAE  conformer  in  the  a-solid.  A  second  crystalline  form  of  RDX  exists  (P-RDX),  but  it  is 
extremely  unstable,  and  no  direct  experimental  structural  information  is  available  for  this  form  [9]. 
Karpowicz  and  Brill  [9]  recorded  the  IR  spectra  of  both  a-  and  P-RDX,  as  well  as  RDX  in  the  vapor 
phase.  The  IR  spectra  of  a-  and  P-RDX  have  distinct  differences,  and  the  fewer  modes  in  the  P-solid 
suggest  a  higher  molecular  symmetry  than  that  of  a-RDX.  Additionally,  the  P-solid  spectrum  has 
features  that  are  similar  to  vapor-phase  RDX  [9].  Karpowicz  and  Brill  concluded  that  the  molecular 
conformation  of  the  RDX  in  the  P-solid  and  in  the  vapor  phase  has  a  molecular  symmetry  of  C3v  and 
suggested  two  possible  structures  with  this  symmetry  [9].  One  structure  has  the  nitro  groups 
occupying  all  axial  positions  (AAA),  and  one  has  the  nitro  groups  occupying  all  pseudo-equatorial 
positions  (EEE).  The  measurements  did  not  provide  sufficient  information  to  distinguish  between 
the  two  possibilities.  Subsequent  electron  diffraction  experiments  indicated  that  the  AAA  RDX 
conformer  is  consistent  with  the  diffraction  patterns,  and  structural  parameters  were  obtained  from 
fitting  models  to  the  experimental  data  [10]. 

We  compare  ab  initio  and  nonlocal  DFT  predictions  of  structural  parameters  and  vibrational 
frequencies  for  the  AAE,  AAA,  and  EEE  conformers  of  RDX  against  the  experimental  information 
[8-10].  Second-order  Moeller-Plesset  (MP2)  [11]  geometry  optimizations  using  the  6-31G*  basis 
set  [12-14]  are  used  to  locate  the  AAE  and  AAA  conformers.  Nonlocal  DFT  geometry 
optimizations  using  the  6-3 1G*  [12-14]  and  6-31 1+G**  [15-16]  basis  sets  and  the  B3LYP  density 
functional  [17-19]  are  performed  for  the  same  conformers  and  compared  against  MP2  to  determine 
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the  magnitude  of  difference  in  theoretical  treatments.  As  shown  hereafter,  the  better  agreement  of 
the  B3LYP/6-31 1+G**  predictions  with  experiment  strongly  suggests  that  this  level  is  sufficient  to 
accurately  determine  stable  structures  for  the  RDX  conformers.  Normal-mode  analyses  are  used  to 
characterize  each  stable  point,  and  IR  vibrational  spectra  are  simulated  for  comparison  with 
experiment.  The  spectra  are  simulated  by  fitting  the  predicted  IR  intensities  to  Lorentzian  functions 
with  bandwidths  arbitrarily  set  to  8.  All  calculations  are  performed  using  the  Gaussian  94  set  of 
quantum  chemistry  programs  [20].  All  geometry  optimizations  meet  or  exceed  the  default 
convergence  criteria  assigned  by  Gaussian  94  [20].  The  DFT  calculations  use  the  default  grid  size 
given  in  Gaussian  94  [20]. 

2.  Results  and  Discussion 

2.1  AAE  Structural  Data.  Table  1  lists  the  geometric  parameters  of  the  AAE  conformer 
predicted  with  various  theoretical  methods  and  provides  a  comparison  with  parameters  obtained 
from  neutron  diffraction  measurements  of  a-RDX  [8].  The  labeling  of  the  atoms  in  the 
two-dimensional  projection  of  the  RDX  molecule  in  Figure  1  is  consistent  with  the  labeling  of  the 
internal  coordinates  in  Table  1.  The  angle  0;  shown  for  the  AAA  conformer  in  Figure  1  is  the  angle 
between  the  C(l)-N(2)-C(6)  plane  and  the  plane  containing  the  C(l)-N(3)-N(10)-C(6)  atoms.  Angle 
02  denotes  the  angle  between  the  planes  containing  the  N(3)-C(7)-N(10)  atoms  and  the  C(l)-N(3)- 
N(10)-C(6)  atoms,  respectively.  The  angle  6  is  the  angle  between  the  plane  of  the  C(l)-N(2)-C(6) 
atoms  and  the  N(2)-N(9)  bond.  Of  the  three  theoretical  treatments,  the  B3LYP/6-311+G** 
predictions  have  the  smallest  overall  deviation  from  experiment  for  the  bonds  and  bond  angles,  and 
the  MP2/6-31G*  has  the  largest  overall  deviation  from  experiment.  The  largest  differences  between 
theoretical  predictions  at  the  three  levels  presented  here,  and  the  experimental  determinations  are  in 
the  N-N  bond  lengths  and  the  C-N-N  angles.  All  of  the  theoretical  predictions  overestimate  the  N-N 
bond  lengths  by  ~2.5-4.5%,  and  the  MP2/6-31G*  calculations  overestimate  the  N-0  bond  lengths 
by  2.1-2.6%.  All  theoretical  methods  also  overestimate  the  C(l)-H(4)  bond  by  ~4%.  The  B3LYP/ 
6-3 1 1+G**  predictions  of  the  remaining  bonds  are  in  closest  overall  agreement  with  experiment. 
All  of  the  theoretical  methods  underestimate  the  C-N(2)-N  bond  angles  by  ~3.5-6.0%.  The 
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Table  1.  Structural  Parameters  of  the  AAE  RDX  Conformer 


Bond 

(A) 

MP2/6-31G* 

B3LYP/6-31G* 

B3LYP/6-3 1 1+G** 

Experiment  [8] 

C(l)-N(2) 

1.4700 

1.4733 

1.4748 

1.464 

C(l)-N(3) 

1.4484 

1.449 

1.4488 

1.443 

C(7)-N(3) 

1.4618 

1.4621 

1.4628 

1.468 

C(7)-N(10) 

1.4618 

1.4621 

1.4628 

1.458 

C(6)-N(10) 

1.4484 

1.449 

1.4488 

1.440 

C(6)-N(2) 

1.4700 

1.4733 

1.4748 

1.450 

C(l)-H(4) 

1.1001 

1.0994 

1.0974 

1.058 

C(l)-H(5) 

1.0846 

1.0841 

1.0827 

1.092 

C(7)-H(17) 

1.0864 

1.0851 

1.0839 

1.085 

C(7)-H(18) 

1.0936 

1.0938 

1.0918 

1.087 

C(6)-H(15) 

1.1001 

1.0994 

1.0974 

1.088 

C(6)-H(16) 

1.0846 

1.0841 

1.0827 

1.075 

N(2)-N(9) 

1.4105 

1.4022 

1.4051 

1.351 

N(3)-N(8) 

1.4370 

1.4317 

1.4335 

1.392 

N(10)-N(19) 

1.4360 

1.4317 

1.4334 

1.398 

N(9)-0(13) 

1.236 

1.2258 

1.2195 

1.209 

N(9)-0(14) 

1.236 

1.2258 

1.2195 

1.233 

N(8)-0(ll) 

1.2319 

1.2202 

1.2137 

1.203 

N(8)-0(12) 

1.2324 

1.2211 

1.2148 

1.207 

N(19)-O(20) 

1.2319 

1.2202 

1.2137 

1.201 

N(19)-0(21) 

1.2324 

1.2211 

1.2148 

1.205 

Angle 

(°) 

MP2/6-31G* 

B3LYP/6-31G* 

B3LYP/6-31 1-t-G** 

Experiment 

6i 

54.14 

50.05 

50.87 

53.31 

02 

42.74 

41.44 

41.44 

43.73 

6 

-27.91 

-24.39 

-24.08 

-12.59 
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Table  1.  Structural  Parameters  of  the  AAE  RDX  Conformer  (continued) 


Angle  MP2/6-31G*  B3LYP/6-31G*  B3LYP/6-311+G**  Experiment 

(°)  

N(2)-C(l)-N(3)  108.84  109.26  109.09  107.8 

N(2)-C(l)-H(4)  111.36  110.90  110.64  109.9 

N(2)-C(l)-H(5)  109.61  109.38  109.57  UOO 

N(3)-C(l)-H(4)  106.90  107.55  107.49  108.0 

N(3)-C(l)-H(5) _ 110.19 _ 109.95 _ 110.30 _ 110.0 

H(4)-C(l)-H(5)  109.89  109.78  109.73  111.0 

N(3)-C(2)-N(10)  113.38  112.71  112.45  111.7 

N(3)-C(2)-H(17)  109.91  109.76  109.99  110.1 

N(3)-C(2)-H(18)  106.80  107.21  107.24  106.9 

N(10)-C(2)-H(17)  109.91  109.76  109.99  110.7 

N(10)-C(2)-H(18)  106.80  107.22  107.24  107.2 

H(17)-C(2)-H(18)  109.93  110.09  109.84  110.1 

N(10)-C(6)-N(2)  108.84  109.26  109.09  108.4 

N(10)-C(6)-H(15)  106.90  107.55  107.49  107.4 

N(10)-C(6)-H(16)  110.19  109.95  110.30  HU 

N(2)-C(6)-H(15)  111.36  110.90  110.64  M6 

N(2)-C(6)-H(16)  109.61  109.38  109.57  HL3 

H(15)-C(6)-H(16)  109.89  109.78  109.73  100 

C(6)-N(2)-C(l)  112.76  115.03  114.52  115.1 

C(6)-N(2)-N(9)  113.65  115.09  115.42  1209 

C(l)-N(2)-N(9)  113.65  115.09  115.42  U9J 

N(2)-N(9)-Q(13)  116.52 _ 116.60  116.64  117.2 

N(2)-N(9)-Q(14)  116.52 _ 116.60 _ 116.64  117.8 

Q(13)-N(9)-Q(14)  126.91  126.78  126.68  125^0 

C(l)-N(3)-C(7)  113.88  115.47  115.52  114.6 
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Table  1.  Structural  Parameters  of  the  AAE  RDX  Conformer  (continued) 


MP2/6-31G* 

B3LYP/6-31G* 

B3LYP/6-31 1+G** 

Experiment 

N(3)-N(8)-0(12) 

116.66 

116.63 

116.65 

116.8 

0(11)-N(8)-0(12) 

127.37 

127.11 

127.11 

125.7 

C(7)-N(10)-C(6) 

113.88 

115.47 

115.52 

114.8 

C(7)-N(10)-N(19) 

114.46 

116.51 

116.91 

117.5 

C(6)-N(10)-N(19) 

113.83 

116.05 

116.43 

115.6 

N(  1 0)-N(l  9)-O(20) 

115.78 

116.11 

116.10 

117.3 

N(10)-N(19)-O(21) 

116.66 

116.63 

116.65 

117.0 

O(20)-N(  1 9)-0(2 1 ) 

127.37 

127.11 

127.11 

125.5 

MP2/6-31-G*  predictions  of  the  remaining  C-N-N  bonds  underestimate  the  values  by  1. 5-2.8%, 
while  the  B3LYP  predictions  are  within  1%  of  experiment.  The  predicted  angles  0,  and  02,  which 
are  indicative  of  the  deviation  of  the  ring  atoms  from  planarity,  are  within  3.1  and  2.4%  or  less, 
respectively,  from  the  experimental  values.  The  theoretical  predictions  of  6  (which  measures  the  tilt 
of  the  N-N  bond  away  from  the  CNC  plane)  disagree  with  experiment  by  1 1.5°  to  15.3°  and  could 
be  due  to  crystal-field  effects.  The  B3LYP/6-31 1+G**  prediction  of  6  is  in  closest  agreement  with 
experiment.  The  overall  good  agreement  in  the  geometries  is  remarkable,  since  the  neutron 
diffraction  information  was  determined  from  molecules  in  the  crystal  state,  and  the  theoretical 
calculations  involve  a  single  RDX  molecule.  The  crystal  field  does  not  significantly  distort  the 
molecule  from  Cs  symmetry. 

2.2  AAA  and  EEE  Structural  Data.  The  AAA  and  EEE  conformers  suggested  as  possible 
structures  for  RDX  in  the  vapor  and  P-solid  phases  both  have  C3v  symmetry  [9],  Theoretical 
predictions  of  geometric  parameters  for  both  conformers  are  given  in  Table  2,  along  with  the 
experimental  information  for  vapor-phase  RDX.  In  this  table,  the  individual  geometric  parameters 
are  given  along  with  the  averages  of  symmetry  equivalent  parameters.  The  averages  are  compared 
against  the  experimental  numbers.  The  theoretical  predictions  of  the  bond  lengths  for  both 
conformers  at  all  levels  are  within  1%  of  the  experimental  result.  The  agreement  of  bond  angles  with 
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Table  2.  Structural  Parameters  of  the  AAA  and  EEE  RDX  Conformers 


= 

AAA 

EEE 

Expt. 

MP2 

B3LYP 

B3LYP 

6-3 1G* 

6-3 1G* 

6-31 1+G** 

6-3 1G* 

1 6-31 1+G** 

CN 

C(l)-N(2) 

1.4586 

1.4603 

1.4607 

1.4612 

1.4612 

C(6)-N(2) 

1.4586 

1.4601 

1.4608 

1.4606 

1.4612 

C(6)-N(10) 

1.4586 

1.4604 

1.4607 

1.4596 

1.4604 

C(7)-N(10) 

1.4586 

1.4599 

1.4605 

1.4604 

1.4605 

C(7)-N(3) 

1.4586 

1.4606 

1.4607 

1.4575 

1.4593 

C(l)-N(3) 

1.4586 

1.4606 

1.4605 

1.4577 

1.4598 

<CN> 

1.4586 

1.4603 

1.4607 

1.4595 

1.4604 

1.464 

NN 

N(2)-N(9) 

1.4280 

1.4230 

1.4237 

1.4011 

1.4043 

N(10)-N(19) 

1.4280 

1.4223 

1.4236 

1.4010 

1.4041 

N(3)-N(8) 

1.4280 

1.4223 

1.4229 

1.3987 

1.4030 

<NN> 

1.4280 

1.4225 

1.4234 

1.4003 

1.4038 

1.413  | 

II_ NO 

N(8)-0(ll) 

1.2331 

1.2219 

1.2157 

1.2254 

1.2189 

N(8)-0(12) 

1.2331 

1.2221 

1.2156 

1.2255 

1.2189 

N(9)-0(13) 

1.2331 

1.2219 

1.2155 

1.2251 

1.2187 

N(9)-0(14) 

1.2331 

1.2218 

1.2156 

1.2252 

1.2188 

N(19)-O(20) 

1.2331 

1.2219 

1.2157 

1.2252 

1.2188 

N(19)-0(21) 

1.2331 

1.2222 

1.2156 

1.2252 

1.2187 

<NO> 

1.2331 

1.2220 

1.2156 

1.2253 

1.2188 

1.213 

Note:  Bond  lengths  in  angstroms  and  angles  in  degrees. 
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Table  2.  Structural  Parameters  of  the  AAA  and  EEE  RDX  Conformers  (continued) 


N(2)-C(6)-N(10) 


N(2)-C(l)-N(3) 


N(3)-C(7)-N(10) 


<NCN> 


MP2 


6-3 1G* 


AAA 


B3LYP 


6-3 1G*  6-31 1+G** 


CH 


EEE 


B3LYP 


6-3 1G*  6-31 1+G** 


C(l)-H(4) 

1.0947 

1.0947 

1.0929 

1.1041 

1.1022 

C(6)-H(15) 

1.0947 

1.0947 

1.0928 

1.1043 

1.1023 

C(7)-H(18) 

1.0947 

1.0948 

1.0928 

1.1043 

1.1022 

C(l)-H(5) 

1.0866 

1.0854 

1.0843 

1.0846 

1.0830 

C(6)-H(16) 

1.0866 

1.0854 

1.0843 

1.0844 

1.0828 

C(7)-H(17) 

1.0857 

1.0856 

1.0843 

1.0845 

1.0829 

<CH> 

1.0905 

1.0901 

1.0886 

1.0944 

1.09269 

NCN 


113.64 

112.71 

112.47 

106.07 

105.92 

113.64 

112.79 

112.43 

105.97 

105.82 

113.64 

112.72 

112.46 

105.92 

105.73 

113.64 

112.74 

112.45 

105.99 

105.82 

C(l)-N(2)-C(6)  114.08  115.50 


C(6)-N(10)-C(7)  114.08  115.67 


C(l)-N(3)-C(7)  114.08  115.62 


<CNC>  1114.08  115.60 


C(l)-N(2)-N(9) 


C(l)-N(3)-N(8) 


C(6)-N(2)-N(9) 


C(6)-N(10)-N(19) 


116.08 


116.07 


116.08 


116.08 


117.45 


117.43 


117.40 


117.46 


CNC 


115.68 


115.72 


115.80 


115.73 


CNN 


117.94 


118.10 


117.94 


117.99 


117.62 


117.44 


117.88 


117.65 


115.13 


115.89 


115.09 


115.31 


116.95 


116.90 


117.12 


116.99 


115.56 


116.01 


115.53 


115.69 
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Table  2.  Structural  Parameters  of  the  AAA  and  EEE  RDX  Conformers  (continued) 


AAA 

EEE 

Expt. 

MP2 

B3LYP 

B3LYP 

6-3 1G* 

6-3 1G* 

6-31 1+G** 

6-3 1G* 

6-31 1+G** 

C(7)-N(10)-N(19) 

116.08 

117.51 

118.00 

115.31 

115.72 

C(7)-N(3)-N(8) 

116.08 

117.48 

118.10 

115.93 

116.01 

<CNN> 

116.08 

117.46 

118.01 

115.44 

115.75 

116.3 

1  ONO 

0(11)-N(8)-0(12) 

127.24 

127.00 

126.97 

126.95 

126.83 

0(13)-N(9)-0(14) 

127.24 

127.04 

126.98 

126.98 

126.84 

O(20)-N(19)-O(21) 

127.24 

127.00 

126.99 

126.98 

126.85 

<ONO> 

127.24 

127.01 

126.98 

126.97 

126.84 

125.5 

|  HCH 

H(4)-C(l)-H(5) 

109.89 

110.09 

109.79 

109.39 

109.56 

H(17)-C(7)-H(18) 

109.89 

110.02 

109.81 

109.41 

109.59 

H(15)-C(6)-H(16) 

109.89 

110.09 

109.79 

109.37 

109.55 

<HCH> 

109.89 

110.07 

109.80 

109.39 

109.57 

105.1 

$ 

-0.03 

-0.01 

-0.07 

-0.01 

0.07 

19.1 

0i 

42.25 

41.75 

42.09 

51.36 

52.45 

33.9 

02 

42.07 

40.48 

40.65 

45.80 

47.00 

6 

23.37 

19.70 

18.41 

-21.01 

-21.38 

19.9 

Y 

346.24 

350.50 

351.75 

348.54 

348.56 

356.3 

Note:  Bond  lengths  in  angstroms  and  angles  in  degrees. 


experiment,  however,  is  not  as  good  for  the  C3v  conformers  as  it  was  for  the  AAE  structure.  The 
largest  disagreement  between  the  calculated  and  experimental  values  is  in  the  CNC  angles.  The 
B3LYP/  6-31 1+G**  predictions  for  AAA  and  EEE  are  115.7°  and  117.0°,  respectively,  which 
underestimate  the  experimental  value  of  123.7°  by  8.0°  to  6.7°.  The  agreement  of  the  MP2/6-31G* 
AAA  prediction  of  this  angle  with  experiment  is  even  worse  (~10°).  It  is  notable  that  the  theoretical 
predictions  of  the  CNC  angle  for  these  conformers  are  closer  to  the  experimental  value  for  the  AAE 
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conformer  (-1 14.8°).  Shiskov  et  al.  [10]  note  a  large  change  in  the  CNC  angles  between  the  gas 
phase  and  a-solid  AAE  conformers,  but  attribute  the  difference  to  crystal-field  effects  for  the  AAE 
conformer.  As  previously  mentioned,  this  is  not  supported  by  the  current  results  where  theory  (for 
the  isolated  AAE  molecule)  predicts  an  AAE  structure  in  close  agreement  with  the  experimental 
crystal  structure.  In  addition,  the  theoretical  predictions  for  the  CNC  angles  are  quite  similar  for  all 
three  conformers  and  agree  with  the  experimental  AAE  crystal  CNC  angle  to  within  2.8°  or  less. 
Theoretical  predictions  of  the  angle  0X  are  greater  than  experiment  by  ~8°  and  ~18°  for  the  AAA 
and  EEE  conformers,  respectively.  These  indicate  that  experiment  predicts  ring  structures  closer  to 
planarity  than  either  of  the  theoretical  C3v  structures,  although  the  AAA  values  are  closer  to 
experiment  than  the  EEE.  We  have  also  included  the  predicted  values  for  y,  which  are  included  in 
the  experimental  paper  and  defined  as  the  sum  of  the  three  bond  angles  involving  the  ring  nitrogen 
atom.  This  value  reflects  the  degree  to  which  the  ring  nitrogen  is  coplanar  with  its  three  attached 
neighbors.  This  parameter  for  both  conformers  at  all  levels  is  in  close  agreement  with  experiment. 

The  geometric  parameter  <j),  defined  by  Shishkov  et  al.  [10]  as  the  torsional  angle  about  the  N-N 
bond,  was  reported  to  be  19.1  °.  The  value  for  the  optimized  structures  calculated  at  all  levels  in  this 
study  is  0°,  indicating  a  geometry  in  which  the  “C. .  .C  and  O. .  .0  lines  of  the  C2N-N02  fragment 
are  coplanar”  [6].  A  geometry  optimization  of  the  AAA  conformer  at  the  B3LYP/6-31G*  level  was 
attempted  in  which  the  starting  geometry  had  the  angle  <|)  set  to  19.3°.  The  value  of  the  angle  4>  in 
the  resulting  optimized  structure  was  0°,  and  the  remaining  geometric  parameters  are  equal  to  those 
given  in  Table  2.  The  discrepancy  between  theory  and  experiment  could  be  related  to  the  fitting 
procedure  used  in  the  experimental  analysis.  In  the  experimental  analysis,  structural  models  were 
assumed,  each  of  which  was  described  by  a  set  of  variables  that  would  be  parameterized  to  provide 
best  agreement  with  the  electron  diffraction  measurements.  Three  assumptions  were  made  about 
structural  relationships  that  were  held  fixed  throughout  the  fitting  procedures:  the  CH2  moieties  have 
local  C2v  symmetry,  the  CNN  angles  within  a  molecule  are  equal,  and  the  N02  geometry  is  planar. 
The  final  geometry  was  obtained  by  an  iterative  fitting  procedure  in  which  one  parameter  was 
optimized  at  a  time  when  the  remaining  parameters  were  held  fixed.  Initially,  all  parameters  were 
assigned  starting  values  taken  from  the  literature  for  similar  compounds.  After  the  iterative 
parameter  refinement,  a  final  simultaneous  least-squares  fit  of  the  parameter  set  was  performed, 
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resulting  in  the  reported  values.  Such  an  iterative  procedure  could  be  subject  to  convergence  to  a 
local  minimum  in  the  parameter  space.  This  is  a  possible  explanation  for  the  difference  between  the 
theoretical  and  experimental  value  of  <J>. 

The  information  gained  by  the  bond  lengths  and  angles  does  not  provide  sufficient  information 
to  distinguish  between  the  C3v  conformers,  nor  do  the  angles  and  <J>2.  However,  the  predicted 
values  of  6  for  AAA  are  almost  in  exact  agreement  with  the  experimental  value,  while  the  EEE 
predictions  are  off  by  41  Therefore,  based  primarily  on  6  and  less  so  on  0lf  the  theoretical 
calculations  clearly  support  the  conclusions  obtained  from  the  experiment  [10];  namely,  that  the 
vapor-phase  structure  of  RDX  is  consistent  with  the  nitro-group  arrangement  in  the  AAA  conformer. 

2.3  Vibrational  Spectra.  Harmonic  vibrational  frequencies  for  the  three  conformers  were 
determined  through  normal-mode  analyses;  each  conformer  had  six  0  frequencies,  and  the  remaining 
frequencies  were  real,  in  contrast  to  the  SCF/4-21G  results  [7].  Table  3  provides  the  calculated 
harmonic  vibrational  frequencies  at  die  MP2/6-31G*  andB3LYP/6-311+G**  levels,  corresponding 
IR  intensities  (in  esu2-cm2)  and  symmetry  assignments  for  each  mode  for  comparison  with  the 
experimental  assignments.  Assignments  identifying  the  nature  of  the  vibrational  modes  are  given 
in  Table  4.  The  symmetry  assignments  correspond  to  the  B3LYP/  6-31 1+G**  results,  which  appear 
to  give  the  best  reproduction  of  the  experimental  spectra,  as  shown  hereafter.  It  is  hoped  that  these 
assignments  will  be  of  help  to  experimentalists  in  interpreting  observed  spectra. 

Simulated  spectra  based  upon  B3LYP/6-3 1 1 +G*  *  IR  intensities  and  frequencies  for  the  three 
conformers  are  compared  against  experimental  IR  spectra  in  Figure  2.  The  vibrational  frequencies 
in  this  figure  have  been  reduced  by  3%. 

It  is  clear  that  the  simulated  IR  spectrum  for  the  AAA  conformer  has  several  features  that  are 
similar  to  the  experimental  vapor  and  P-solid  phase  IR  spectra  in  the  mid-infrared  (IR)  region, 
particularly  between  1,100  cm'1  and  1,650  cm'1  [9].  The  theoretical  IR  spectrum  between  1,100  and 
1,500  cm'1  for  the  EEE  conformer  shows  a  band  pattern  in  much  poorer  agreement  with  the 
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Table  3.  Theoretical  and  Experimental  Vibrational  Frequencies  of  RDX  Conformers 


EEE  Conformer  1! 

Theoretical  EEE 

* 

* 

o 

+ 

1-H 

2 
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co 

« 

w 

W 

W 

ffl 

i 

w 

W 

W 

b 

1 

h 

E 

CN 

cs 

t- 

^H 

E 

NO 

CS 

1-H 

VO 

CS 

1-H 

E 

E 

OS 

i-H 

OS 

^H 

0 

¥ 

0 

$ 

to 

00 

3,203 

3,200 

3,199 

s 

os 

cs 

2,959 

2,958 

CO 

VO 

1-H 

r-* 

CO 

vq 

1—H 

1,619 

1,522 

1,510 

8 

n 

1-H 

1,431 

0 

9 

i-H 

8 

CO 

1-H 

1,350 

vo 

CO 

cq 

i-H 

AAA  Conformer 

i 

c a 

i 

CQ. 

Ref.  [9] 

IR 

3,075  w 

1 

LJ 

Ed 

1,588  s,sh 

1,458  w 

1,441  m,b 

1,419  m 

£ 

CO 

00 

CO 

1 

i 

Vapor 

Ref.  [9] 

IR 

3,065  vw 

1 

1 

1,584  vs 

1 

1,444  m 

1,420  m 

1,374  m 

1 

1 

Theoretical  AAA 

B3LYP/6-31 1+G** 

hr.  Rep. 

w 

W 

< 

w 

w 

w 

< 

B 

w 

w 

W 

W 

W 

w 

E 

Int. 

E 

00 

1-H 

0 

wo 

1-H 

- 

1,123 

ZZVl 

O 

220 

VO 

00 

00 

Os 

to 

vo 

to 

CO 

0 

^H 

1 

i 

3,194 

3,192 

3,192 

3,070 

3,064 

3,064 

00 

to 

VO 

1-H 

1,658 

1,627 

1,482 

1,466 

1,464 

1,403 

1,402 

N- 

00 

CO 

1,381 

CO 

vo 

cq 

MP2/6-31 

Int 

1-H 

00 

CS 

O 

O 

00 

vo 

^H 

00 

to 

O 

OS 

00 

VO 

CN 

CS 

00 

00 

0 

O 

O 

I 

3,265 

3,265 

3,265 

3,144 

3,138 

3,138 

VO 

CO 

00 

VO 

CO 

00^ 

1-H 

1,828 

1,530 

1,508 

00 

0 

vq 

hH 

1,420 

1,420 

1,407 

1,407 

1,373 

O 

,5 

a-solid 

Ref.  [21] 

Raman 

3,075  m 

1 

Ed 

GA 

1-H 

8 

CO 

6 

1 

CS 

1 

£ 

os 

to 

£ 

0 

r^ 

vq 

i-h 

* 

00 

CO 

vq 

1-H 

1,508  vw 

1,456  w 

* 

<r> 

S? 

1-H 

iC 

CA 

CS 

S' 

1,387  w 

1,377  w 

* 

co 

i-H 

1,309  s 

g 

CA 

to 

b 

CO 

1 

CA 

co" 

6 

1-H 

8 

co" 

p 

1 

CA 

> 

00 

os 

iq 

r-H 

CA 

> 

CO 

r- 

vq 

1-H 

CA 

> 

8 

vq 

CA 

CS 

CO 

vq 

1-H 

CA 

OS 

3 

i-H 

E 

it 

3 

r-H 

E 

<r> 

1-H 

CA 

Os 

00 

CO 

1,377 

T 

cs 

to 

cq 

1-H 

1,310  s 

Ref.  [9] 

Ed 

1 

8. 

Ed 

1 

1 

1,592  s 

1,576  s 

1,539  m 

s 

CO 

CO 

vq 

e 

00 

3 

^H 

E 

SO 

E 

co 

3 

1-H 

6 

00 

00 

CO 

1-H 

1 

s 

to 

cq 

^H 

E 

cq 

i-H 

I 

•a 

O 

§ 

* 

* 

? 

1-H 

1-H 

Q 

(2 

J 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

'< 

< 

< 

< 

< 

'< 

< 

I 

C 

o 

1-H 

00 

T-H 

VO 

to 

r-* 

1-H 

1-H 

VO 

s 

to 

to 

CO 

to 

1— H 

CO 

CS 

1-H 

¥ 

CN 

Os 

1-H 

VO 

00 

¥ 

i-H 

2 

CO 

PQ 

1 

3,206 

3,205 

3,199 

3,081 

3,016 

3,015 

1,668 

vq 

CO 

8 

vo 

§ 

i-H 

1,480 

00 

O 

3 

1— H 

1,406 

1,374 

1,362 

CS 

VO 

cq 

O 

§ 

o 

1-H 

CO 

1 

VO 

|  Int. 

00 

3 

CO 

1-H 

^H 

CS 

1— H 

Os 

VO 

r- 

cs 

r- 

cs 

1-H 

1-H 

vo 

1-H 

1-H 

124 

3 

1H 

129 

0 

| 

3,277 

3,277 

3,272 

3,155 

3,083 

3,083 

00 

1,837 

1,816 

1,547 

1,529 

00 

1-H 

vq 

i-H 

1,445 

1,420 

1,395 

1,382 

1,376 1 

_ 

ss 

□ 

1 

Bl 

Ql 

ESI 

□I 

Bl 

Bl 

Bj 

Bj 

Bj 

Bl 

D 

Bl 

Bl 

Bl 

B 

Bl 

m 
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Notes:  vs  -  very  strong,  s  -  strong,  m  -  medium,  b  -  broad,  sh  -  shoulder. 

Frequencies  and  corresponding  IR  intensities  are  in  cm'1  and  esu2-cm2,  respectively. 
Symmetry  assignments  correspond  to  the  B3LYP/6-3 1 1+G**  results  only. 


Table  3.  Theoretical  and  Experimental  Vibrational  Frequencies  of  RDX  Conformers  (continued) 


EEE  Conformer 

Theoretical  EEE 

B3LYP/6-31 1+G**| 

I 

w 

w 

pq 

pq 

pq 

pq 

pq 

pq 

< 

pq 

w 

pq 

Int. 

T-H 

ON 

© 

ON 

VO 

CN 

vo 

t-H 

CO 

VO 

T-H 

CN 

t-H 

© 

VO 

vo 

t-H 

Tf 

00 

CO 

$ 

t-H 

209 
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wo 

CN 

t 

s 

1,334 

& 

CN 

t*H 

£ 

CN 

t»»H 

Tf 

CN 

t*H 

tC 

CN 

t*H 
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o 

00 

^H 

§ 

ON 

© 

t-H 
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© 

T-H 

On 

00 

00 

00 

00 

OO 

00 

819 
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P-solid 

Ref.  [9] 

g 
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1 
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g 
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1 
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1 
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1 
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r- 
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Notes:  vs  -  very  strong,  s  -  strong,  m  -  medium,  b  -  broad,  sh  -  shoulder. 

Frequencies  and  corresponding  IR  intensities  are  in  cm"1  and  esu2-cm\  respectively. 
Symmetry  assignments  correspond  to  the  B3LYP/6-31 1+G**  results  only. 


Table  3.  Theoretical  and  Experimental  Vibrational  Frequencies  of  RDX  Conformers  (continued) 
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Notes:  vs  -  very  strong,  s  -  strong,  m  -  medium,  b  -  broad,  sh  -  shoulder. 

Frequencies  and  corresponding  IR  intensities  are  in  cm'1  and  esu2-cm2,  respectively. 
Symmetry  assignments  correspond  to  the  B3LYP/6-31 1+G**  results  only. 


Table  3.  Theoretical  and  Experimental  Vibrational  Frequencies  of  RDX  Conformers  (continued) 


Table  4.  Theoretical  Vibrational  Frequencies  (cm-1)  and  Assignments  of  RDX  Conformers 


<3 

i 

W 

S 

CH  st  (eq)  1 

CH  st  (eq) 

CH  st  (eq) 

CH  st  (ax) 

CH  st  (ax) 

CH  st  (ax) 

to 

o 

£ 

O 

0S 

o 

o 

to 

o 

1 

o 

*3 

CO 

g 

*3 

« 

a? 

u 

•  H 

o 

CO 

g 

3> 

* 

5" 

Sbm.  zhd 

sp 

* 

if 

to 

o 

! 

+ 

4-4 

CO 

1 

4-4 

CO 

+ 

B 

if 

I 

IX, 

o 

04 

r* 

CO 

1. 

CO 

ON 

ON 

CO 

NO 

04 

ON 

£ 

<s 

00 

ON 

04* 

0> 

CO 

NO 

0- 

CO 

so 

•* 

Os 

$ 

t-H 

8 

O 

*n 

8 

NO 

H 

9 

O 

5? 

so 

SO 

CO 

•* 

i-H 

O 

NO 

CO 

* 

SO 

CO 

CO 

* 

AAA  Conformer 

4-» 

I 

■  -4 

3 

? 

S— / 

4H 

CO 

X 

u 

w 

'w' 

4—* 

CO 

5 

? 

N—" 

4-> 

CO 

g 

s 

£ 

'■w' 

4-4 

CO 

g 

? 

4-* 

CO 

g 

? 

'w' 

4-» 

CO 

g 

to 

o 

Z 

o 

4-4 

CO 

O 

A 

£ 

a 

4-4 

CO 

o 

A 

s 

o 

•a 

a 

g* 

* 

a 

g 

-s 

a 

g 

8P 

£ 

if 

& 

* 

g 

B 

g 

B 

g 

SP 

£ 

g 

i 

Tt 

ON 

CO 

Os 

co* 

§ 

co* 

1 

CO 

I 

CO* 

i 

CO 

00 

in 

SO 

00 

NO 

SO 

t-H 

C4 

NO 

04 

so 

$ 

t-H 

1 

t-H 

1 

t-H 

§ 

Tf- 

00 

CO 

t-H 

00 

CO 

t-H 

CO 

so 

CO 

t-H 

AAE  Conformer 

i 

*3 

< 

? 

to 

EC 

0 

f 

'w' 

4-4 

CO 

5 

4-» 

CO 

1 

4-* 

CO 

1 

? 

w 

4-4 

CO 

g 

? 

to 

g 

CO 

o 

£ 

o 

1 

4-4 

CO 

o 

o 

t 

4-4 

CO 

O 

A 

* 

o 

*3 

CO 

tf 

•  t-< 

o 

CO 

if 

o 

CO 

£ 

8P 

& 

g 

gp 

* 

g 

B 

g 

8P 

£ 

if 

4-4 

CO 

§ 

B 

g 

1 

i 

s 

cs 

CO 

O 

04 

CO*4 

ON 

t-H 

CO* 

00 

© 

CO 

so 

o 

CO* 

ID 

t-H 

o 

CO* 

00 

SO 

SO 

SO 

CO 

04 

so 

t-H 

| 

| 

00 

5? 

| 

t-H 

r- 

CO 

t-H 

04 

so 

CO 

t-H 

04 

SO 

CO 

■ 

§ 

t-H 

04 

CO 

in 

so 

0- 

00 

ON 

o 

04 

t-H 

CO 

t-H 

NO 

SO 

0- 

17 


Notes:  Evaluated  at  the  B3LYP/6-31 1+G*  level. 

st  =  stretch,  b  =  bend,  tw  =  twist,  umb  =  umbrella,  rot  =  rotation,  sci  =  scissor. 
The  most  prominent  contribution  to  each  mode  is  listed  first 


Table  4.  Theoretical  Vibrational  Frequencies  (cm1)  and  Assignments  of  RDX  Conformers  (continued) 
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Notes:  Evaluated  at  the  B3LYP/6-3 1 1+G*  level. 

st  =  stretch,  b  =  bend,  tw  =  twist,  umb  -  umbrella,  rot  =  rotation,  sci  =  scissor. 
The  most  prominent  contribution  to  each  mode  is  listed  first 


Table  4.  Theoretical  Vibrational  Frequencies  (cm  *)  and  Assignments  of  RDX  Conformers  (continued) 


Notes:  Evaluated  at  the  B 3LYP/6-3 1 1 +G*  level. 

st  =  stretch,  b  =  bend,  tw  =  twist,  umb  =  umbrella,  rot  =  rotation,  sci  =  scissor. 
The  most  prominent  contribution  to  each  mode  is  listed  first 


Table  4.  Theoretical  Vibrational  Frequencies  (cm1)  and  Assignments  of  RDX  Conformers  (continued) 
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Figure  2.  Simulated  IR  Spectra  at  the  B3LYP/6-311+G**  Level  of  the  AAA,  EEE,  and  AAE 
Conformers.  The  Vibrational  Frequencies  Used  to  Generate  These  Spectra  Are 
Reduced  by  3%.  Experimental  Spectra  From  Karpowicz  and  Brill  [9]  for  a-RDX, 
P-RDX,  and  Vapor-Phase  RDX  Are  Given  for  Comparison. 
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experimental  spectra.  These  spectra  offer  strong  support  for  the  conclusions  of  the  electron 
diffraction  study  that  the  AAA  conformer  is  the  most  probable  structure  of  RDX  in  the  vapor  phase 
[10].  These  spectra  also  show  that  theoretical  treatments  are  useful  in  differentiating  between  two 
possible  conformers  in  the  absence  of  diffraction  data.  We  have  also  provided  a  comparison  of  our 
simulated  AAE IR  spectrum  (frequencies  reduced  by  3%)  with  the  experimental  spectrum  for  the 
a-solid  form  in  Figure  2.  With  the  exception  of  the  bands  around  1,040  cm'1  in  the  experimental 
spectrum,  most  of  the  remaining  experimental  bands  can  be  clearly  assigned  to  a  theoretical 
counterpart.  This  agreement  between  theory  and  experiment  seems  remarkably  good,  considering 
that  the  experimental  spectrum  includes  effects  of  the  crystal  field,  as  well  as  overtones  and 
combination  bands. 

Figure  3  provides  a  comparison  of  simulated  IR  spectra  (using  unsealed  frequencies)  for  the  three 
conformers  at  the  B3LYP  level  using  the  6-31G*  and  6-31 1+G**  basis  sets.  Figure  3  also  includes 
spectra  using  MP2/6-3 1G*  results  for  the  AAE  and  AAA  conformers.  It  is  apparent  that  the  features 
in  the  B3LYP  spectra  are  relatively  insensitive  to  basis  sets.  One  of  the  prominent  discrepancies 
between  the  two  basis  sets  can  be  seen  in  all  three  conformers.  This  is  the  shift  in  the  band  reported 
in  the  1,700-1,800  cm'1  range  at  the  double-zeta  (DZ)  level  to  the  lower-energy  1,600-1,700-cm"1 
range  for  the  triple-zeta  (TZ)  basis.  A  second  prominent  feature  is  evident  in  the  AAE  spectra  and 
is  seen  in  the  relative  intensities  of  the  two  bands  at  1,337  and  1,362  cm'1  from  the  TZ  basis  set. 
These  bands  can  mix  since  both  are  of  A  symmetry,  and  the  6-3 1 1+G**  intensity  of  each  band  is  530 
and  172  esu2-cm2,  respectively.  This  ordering  agrees  with  experiment  (see  Figure  2).  The  relative 
locations  of  these  bands  appear  to  be  reversed  at  the  DZ  level.  It  is  apparent  that  the  MP2/6-3 1G* 
spectrum  for  the  AAE  conformer  is  substantially  different  from  both  B3LYP  spectra.  The 
MP2/6-31G*  prediction  of  the  spectrum  of  the  AAA  conformer  also  has  features  that  differ  from  the 
B3LYP  predictions,  although  not  as  pronounced  as  for  the  AAE  comparison. 

Absolute,  relative,  and  zero-point  energies  of  each  conformer  are  given  in  Table  5.  The  B3LYP 
predictions  using  both  6-3 1G*  and  6-31 1+G**  basis  sets  indicate  that  the  AAE  conformer  has  the 
lowest  energy  of  the  three  conformers,  but  only  by  a  fraction  of  a  kcal/mol.  Similarly,  the 
MP2/6-31G*  predictions  indicate  that  the  AAE  conformer  is  only  slightly  lower  in  energy  than  the 
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Figure  3.  Simulated  IR  Spectra  at  the  B3LYP  Level  Using  the  6-31G*  and  6-311+G**  Basis 
Sets  for  the  AAA,  EEE,  and  AAE  Conformers.  MP2/6-31G*  Spectra  Are  Also 
Included.  Unsealed  Vibrational  Frequencies  Were  Used  to  Generate  These  Spectra. 
DZ  Denotes  the  6-316*  Basis  Set,  and  TZ  Denotes  the  6-311+G**  Basis  Set 
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Table  5.  Absolute  and  Relative  Energies  of  RDX  Conformers 
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AAA  conformer  (by  0.13  kcal/mol).  The  zero-point-corrected  B3LYP/  6-31G*  and  B3LYP/ 
6-31 1+G**  energies  of  the  AAA  conformer  are  0.18  and  0.64  kcal/mol,  respectively,  relative  to 
AAE.  Within  the  level  of  accuracy  for  the  calculations,  the  AAE  and  AAA  conformers  are  identical 
in  their  stability.  The  B3LYP/6-31G*  and  B3LYP/6-31 1-tG**  zero-point-corrected  energies  of  the 
EEE  conformer  are  4.97  and  4.57  kcal/mol,  respectively,  relative  to  AAE.  The  earlier  SCF/4-21G 
calculations  had  the  energy  ordering  reversed  for  the  AAE  and  AAA  conformers  [7].  In  that  study, 
the  AAA  was  more  stable  by  0.6  kcal/mol  [7].  Also,  the  EEE  was  less  stable  than  AAA  by 
7.2  kcal/mol  [7].  Karpowicz  and  Brill  suggest  that  the  intermolecular  forces  of  neighboring  RDX 
molecules  in  the  a-RDX  crystal  are  responsible  for  “the  energetically  unlikely  positioning  of  the  N02 
groups  that  produces  approximately  Cs  molecular  symmetry”  [9].  The  calculations  presented  here 
suggest  that  the  positioning  of  the  N02  groups  in  the  AAE  conformer  is  not  “energetically  unlikely.” 

3.  Conclusions 


Three  conformers  of  the  large  polyatomic  explosive  RDX  have  been  located,  and  their  structures 
and  vibrational  spectra  characterized  with  nonlocal  DFT  treatments  using  DZ  and  TZ  quality  basis 
sets.  The  three  conformers  have  the  N02  oriented  in  either  AAE,  AAA,  and  F.F.F.  arrangements 
relative  to  the  ring.  The  AAA  structure  is  consistent  with  electron  diffraction  results  of  vapor-phase 
RDX  [10],  and  the  AAE  conformer  is  consistent  with  that  of  the  room-temperature-stable  RDX 
crystal  [8].  Additionally,  the  AAA  and  AAE  conformers  are  calculated  using  MP2  theory  for 
comparison  with  the  DFT  predictions.  Of  the  three  levels  of  theoretical  treatment,  the 
B3LYP/6-311+G**  predictions  of  the  geometry  of  the  AAE  conformer  are  in  closest  overall 
agreement  with  experiment,  and  MP2/6-31G*  predictions  are  in  the  poorest  overall  agreement  with 
experiment.  The  B3LYP/6-311+G**  results  for  the  AAE  conformer  are  in  agreement  with 
experiment  to  within  2%  for  all  bond  lengths  with  the  exception  of  a  single  C-H  bond  (3.7%)  and 
the  N-N  bonds  (2.5-4.0%).  The  B3LYP/6-31 1+G**  predictions  of  bond  angles  are  within  1 .6%  of 
experiment  with  the  exception  of  two  C-N-N  angles  (3.6-4.5%).  In  general,  the  B3LYP/6-3 1 1+G** 
level  produces  the  closest  overall  agreement  with  all  available  experimental  data  for  structures  and 
spectra. 
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Structural  parameters  for  the  AAA  conformer  are  in  closer  agreement  to  experiment  [10]  than 
those  predicted  for  the  EEE  conformer.  Additionally,  simulated  IR  vibrational  spectra  of  the  AAA 
conformer  compare  well  with  experimental  spectra  of  vapor-phase  and  P-solid  RDX,  while  the 
simulated  spectrum  of  the  EEE  conformer  did  not.  The  differences  in  predicted  geometries  and 
vibrational  spectra  between  the  AAA  and  EEE  conformers  support  the  experimental  conclusions  that 
RDX  in  the  vapor-  and  P-solid  phases  have  C3v  symmetry  [9]  and  have  the  nitro  groups  arranged  in 
the  AAA  configuration  [10].  In  addition  to  providing  atomic-level  information  about  a  well-studied 
explosive,  the  results  presented  here  provide  another  indication  that  DFT  methods  can  be  applied 
to  large  polyatomic  molecules  with  a  small  computational  cost  and  reliable  results  for  molecular 
structure,  intramolecular  force  fields,  and  vibrational  spectra. 
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